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We study the propagation of gravitational waves (gw) in a uniformly magnetized plasma at 
arbitrary angles to the magnetic field. No a priori assumptions are made about the temperature, 
and we consider both a plasma at rest and a plasma flowing out at ultra-relativistic velocities. In 
the 3+1 orthonormal tetrad description, we find that all three fundamental low-frequency plasma 
wave modes are excited by the GW. Alfven waves are excited by a x polarized GW, whereas the slow 
and fast magneto-acoustic modes couple to the + polarization. The slow mode, however, doesn't 
interact coherently with the GW. The most relevant wave mode is the fast magneto-acoustic mode 
which in a strongly magnetized plasma has a vanishingly small phase lag with respect to the GW 
allowing for coherent interaction over large length scales. When the background magnetic field is 
almost, but not entirely, parallel to the GW's direction of propagation even the Alfven waves grow 
to first order in the GW amplitude. Finally, we calculate the growth of the magneto-acoustic waves 
and the damping of the GW. 
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I. INTRODUCTION 

The interaction of gravitational waves (gw) with a 
magnetized plasma relies on the non-isotropy of the mag- 
netic field and the fact that both GW and electromagnetic 
waves propagate at the speed of light in vacuum. In vac- 
uum, the GW excites an electromagnetic wave propagat- 
ing in the same direction as the GW at the speed of light 
as was first realized by . Propagation through a perfect 
collisionless fluid does not affect a GW 0, II] , but in an 
astrophysical plasma the electromagnetic field is coupled 
to the matter and, indirectly, the GW interacts with both 
the electromagnetic fields and quantities such as density, 
pressure and currents p|. 

For a summary of previous work on the interaction of 
a GW with plasma waves we refer to the introduction of 
P|, referred to from hereon as Paper I, and references 
therein and 6] for more recent work. 

In this paper we study the full problem of a plane 
fronted, monochromatic GW of either + or x polariza- 
tion propagating obliquely through a magnetized colli- 
sionless plasma. Furthermore, we don't specify whether 
the plasma is Poynting flux or matter dominated. Also, 
we allow for relativistic velocities as we want to apply 
our results to an ultra-relativistic force-free wind or jet. 

The outline of this paper is as follows. We start in 
Sect. [H] with a discussion of the Einstein field equations 
(efe) describing both the background curvature due to 
the static energy content and the dynamical interaction 
of space-time with a time dependent energy-momentum 
density. In the geometric optics limit (Sect. IIIB|I the 
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GW can be treated just like photons traveling on null 
geodesies. In Sect.s III Dl - III El we recapitulate how the 
proper reference frame of an observer (the 3 + 1 or- 
thonormal tetrad frame) is defined by taking snapshots 
of four-dimensional space-time and using ordinary vec- 
tor calculus on these spatial hypersurfaces. Also, the 
covariant derivatives and connection coefficients for such 
a non-coordinate frame are summarized. A closed set 
of linearized general relativistic magneto-hydrodynamic 
(mhd) equations is derived in Sect. Hill and solved alge- 
braically in Sect. IIVI We find that all the fundamental 
plasma waves are excited: a + polarized GW excites fast 
magneto- acoustic waves, which was expected from previ- 
ous idealized calculations in J?) and Paper I. In this paper 
we extend this result to a realistic fast mode with both 
electromagnetic and gas properties and find that also the 
slow magneto-acoustic mode is excited (Sect. ITv**E)l . A 
completely new result is that a x polarized GW propa- 
gating at an angle to an ambient magnetic field excites 
Alfven waves in the plasma (Sect. lTvD|) . The space-time 
solutions in the comoving frame are derived in Sect. 
and boosted to an ultra-relativistic wind in the the ob- 
server frame in Sect. lVll In the limit where the phase ve- 
locities of the Alfven and fast mode approach the speed 
of light the waves can interact coherently, which results 
in linear growth of the amplitudes. For the Alfven mode 
this also requires that the angle with the magnetic field 
is very small (Sect. IvTIj) . 

Finally, we investigate the back-reaction on the GW in 
Sect. IVIII and find that as the plasma waves grow, the 
GW is damped with a group velocity that decreases lin- 
early with distance. An intuitive interpretation of some 
of these results is presented in Sect. I Villi and we end 
with conclusions in Sect. ILXl 

Throughout this paper Gaussian geometrized units are 
adopted (G = c = 1). Latin indices a...e stand for 
0, 1, 2, 3, and k for spatial components 1, 2, 3. 
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II. EINSTEIN-MAXWELL EQUATIONS 

In Lorentz gauge, the linearized Einstein field equa- 
tions (efe) for weak gravitational waves interacting with 
a magneto- fluid read @: 

G ab ~ -)ph ab = 8ir5T ab , (1) 

where ST ab is the oscillatory part of the energy- 
momentum tensor. If we further specify to the transverse 
traceless (tt) gauge and consider a GW propagating in 
the z direction, the only independent components of h ab 
are h xx = —h yy = h + and h xy = h yx = h x indicating the 
+ and x GW polarizations respectively. The individual 
propagation equations for the two polarizations are: 

Dh+ = -8n(ST xx - ST yy ), (2a) 
Dh x = ~8ir(6T xy + 6T yx ). (2b) 

As an example, a GW propagating perpendicularly to a 
background magnetic field B° excites perturbations B\ 
that produce an oscillating cross term in the stress-energy 
tensor and take the form Q: 

Uh + {z,t)=AB° x B 1 x {z 1 t), Uh y =Q. (3) 

From Maxwell's equations (coupled to the efe via the 
twice contracted Bianchi identities and consequently con- 
servation of energy-momentum Vf,T ab = 0) one finds 
wave equations for the plasma quantities in which the 
GW appears as a source term: 

&Bl(z,t) cx h+(z,t)B° x , (4) 

where □ is some wave propagator. J5J) together with Q 
self-consistently determine the interaction between the 
GW and the plasma waves (as we will work out in more 
detail in Sect. EDI- 

A. Background curvature 



an example, a background magnetic field comparable to 
the surface field of a neutron star (B° ~ 10 s T) curves 
space on a scale of ~ 10 10 m. If this field would remain 
constant with distance, all the GW energy would be con- 
verted to em energy on a length scale of order 1Z. In 
reality, we will be studying a force-free plasma wind in 
which the magnetic field falls off linearly with distance. 
In this case the background curvature also decreases and 
the interaction length scale is much smaller than the ra- 
dius of curvature. 

The opposite limit of a weak primordial magnetic field 
with an extremely large spatial extent was discussed by 
|3 who argued that it is very difficult to keep the inter- 
action coherent on such a scale in a universe that is not 
a perfect vacuum. 

B. Geometric optics 

The gravitational waves are assumed to be of the form: 
/i+.x ~ 'H(z)e %u ^ z ~ t ' with a slowly varying amplitude 
such that ujH(z) » ^W(z) and ± < 1. This is the 
short wavelength, geometric optics limit where Gws be- 
have as rays moving on null geodesies k a k a = that 
experience dispersion, refraction, lensing etc. The GW 
move in an essentially flat Minkowski background, rj ab , 
so the full metric is g ab = r/ ab + h ab with \h ab \ < 1 and 
self-interactions (of order h 2 ) are negligible. 

We will study the propagation of transverse traceless 
gravitational waves in a magnetized plasma at arbitrary 
angle 9 to a background magnetic field B$. Also we 
don't specify the temperature (or equivalently the pres- 
sure) regime except that the temperature equilibration 
time between the different particle species is short as 
compared to other characteristic timescales to comply 
with a hydrodynamic description. The only other non- 
vanishing plasma quantities in the equilibrium state are 
then the energy density fi and pressure p°. The pass- 
ing GW will excite small (first-order) perturbations in all 
plasma quantities, denoted as for instance: ji = fi° + fi . 



The exact (non-linear) Einstein field equations de- 
scribe the curvature of space-time due to the presence 
of matter and energy. This curvature is described by 
the contracted Riemann tensor with an associated char- 
acteristic length scale, 1Z, given by the magnitude of it's 
largest components as: 

R ab = 8tt (r ab - ^9 ab T^j , (5a) 

^ - (5°) 2 . (5b) 

When the interaction with a magnetic field excites elec- 
tromagnetic waves growing linearly with distance, viz 
B 1 oc B°h + kz we can see from our rough estimate (|5b|) 
that the fraction of GW energy that is converted into elec- 
tromagnetic waves is proportional to B 1 /h+ oc z/lZ. As 



C. No coupling to unmagnetized plasma 

The stress-energy tensor for a homogeneous perfect 
fluid in the rest frame of an observer (u a — (1, 0, 0, 0)) is: 

T ab = {p + p)u a u b + pg ab . (6) 

A linearly polarized GW will produce perturbations in p, 
p and u l of order /i+. x - However, the 8T lJ components 
are all higher order except for the trace ST", which is 
purely gauge. Explicitely, in (J2J T xy = T yx = and 

rpxx _ rpyy _ Q 

Therefore a gravitational wave cannot couple to an un- 
magnetized perfect fluid in linearized theory [3j. Only 
relativistically gyrating particles can interact with a x 
polarized GW through non- linear (p+p)viVj terms [lOllllj . 
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D. Space-time split 

To simplify equations, we will use the so-called 3+1 
split of space-time [lj, Il3 • A time-like observer moving 
with 4-velocity u a perceives space as the 3 dimensional 
hypersurface [H| orthogonal to u a , and u a itself as the 
time axis. We can define u a (x b ) at each point in space- 
time as the direction of time, and define space as the 
snapshots of space at constant time. Subsequently, we 
can split equations into their space and time components 
by using the parallel and orthogonal projection operators 
U a b = —u a Ub and H a b = gab + u a Ub, respectively, with 
U%+H% = 6%. 

As an example, the covariant electromagnetic Faraday 
tensor F ab (and it's dual T ab = \e abcd F cd ) can be sim- 
plified by splitting it into its time and space components: 

ab {U a c U b d + H a c H b d )F cd 
= u a E b - E a u b + e abc B r ■■ 



Covariant derivatives are defined as: V a Ti 



br- 



ibe ■ 



F 



bc B r 



(7) 



where we have defined B a = ^e abc F bc and E a = F ab Ub 
that reduce to the magnetic and electric field, respec- 
tively, in the rest frame of an observer. The last equal- 
ity in @ is only valid in the ideal mhd approximation 
E a = 0. Since E a = e abc u b B c or E = -v x B the require- 
ment that the comoving electric field vanishes replaces 
Ohm's law. Similarly, the covariant Maxwell equations 
can be split into their s pac e and time components. For 
a derivation we refer to 



E. Proper reference frame 

In describing the interaction of a GW with a plasma, 
one has two choices for the reference frame. One is the 
transverse-traceless coordinate frame, discussed in the 
previous sections, which is tuned to a GW with metric: 



9tt 





+ h+(z,t) 
h x (z,t) 1 






h x {z,t) 

-h + (z,t) 

1 



The natural reference frame of an observer, however, is 
an orthonormal tetrad frame (onf). The basis vectors 
that remain orthogonal in the presence of a tt plane 
polarized GW are: 



e - (|;,0,0,0), 

e i = (0>(l-¥)&.-¥&.°)> 



(o, -¥£.(1 + ^)4.0), 

(0,0,0, 



(9) 



e 2 
e 3 



a; J' 



where the partial derivatives reflect the notion that in 
curved space-time, tangent vectors and partial deriva- 
tives are equivalent. With respect to these basis vectors 
the metric is r/ ab . 



T d ba T dc — T d ca T bd , where the connection coefficients are 
linear combinations of the commutation functions "/ a bc 
and not derivatives of the metric as in a coordinate frame 

EM. 



ribdl d ca 



+ Vcdl d ab), (10) 



and are skew in the first two indices, r7 M c = 0. 

There seems to be some inconsistency in the literature 
on the explicit form of the connection coefficients. In a 
coordinate frame, the connection coefficients are given by 
the Christoffel symbols. For the metric JSJ the Christof- 
fel symbols are oc h/2 and have 12 non- vanishing com- 
ponents for each polarization. The explicit form of these 
components are given in and 0, but it should be 
noted that in [14( one of course also has the symmetric 
components, and in Q r° 10 , r° 01 should be r 1 10 , 

In the orthonormal tetrad frame the connection coef- 
ficients only have 8 non-vanishing components for each 
polarization: 



[02] 2 



[01] 1 



T[20]l — — r[01]2 

where T 



l dh+ 
2~df 

1 dh x 
2~ 



[32] 2 



-r [3 i] 



ldh 
2 



dz 



[13]2 



= -r 



[32] 1 



1 dh 
2 



dz 



Of 

[oiji stands for Ton = — Tioi etc. 
The Einstein field equations and the equations of 
geodesic deviation are derived from the Riemann cur- 
vature tensor, that in the onf is given to first order in 
h+,x by: 



R a bcd = e L T\ d -e d T\ c + 0[h 2 



(11) 



The Ricci tensor R ab is just a contraction of (|llfl and to 
first order reduces to the same form as in the tt coordi- 
nate frame. This means that the Einstein field equations 
in the proper reference frame also have the same form as 
in (we will use this in Sect. IVlH when discussing the 
damping of the Gw). 

The driving force of a GW on a test particle is also 
described through the Riemann tensor in the form of the 
equations of geodesic deviation: 



dV 



-R i0 jox j 



h + 


hx 














■ 


. (12 








o / 







where h + 



: J^/i-|_ )X . This equation will be important 
for our interpretation of the interaction with a magnetic 
field in Sect. 171111 



III. GENERAL RELATIVISTIC MHD 

Throughout this paper we assume that the MHD ap- 
proximation is valid, viz that the plasma is a collision- 
less one-fluid with negligible viscosity, resistivity and heat 
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flow. In contrast to some definitions of the mhd approx- 
imation, however, we do allow for relativistic velocities 
which means that displacement currents cannot be ne- 
glected and a generalized definition of the Alfven velocity 
is needed fSect. HV}) . We don't restrict the angle between 
the GW propagation and the background magnetic field, 
but without loss of generality choose it to lie in the x-z 
plane. 



A. Coupling to the electromagnetic field 

A gravitational wave propagating through a uniform 
magnetic field B° produces a Lorenz force F£ given by 
the projection of Ampere's law, V b F ah — j°, perpendic- 
ular to B°: 

i -i „o (S°-V)B 1 / B° ■ B 1 
Ft = j 1 x B = - ^ V ' 



4tt 



4tt 



|B°| 2 dv 1 B° d i j E xB° 

Faraday's law, V b T ah = 0, is given in the 3 + 1 split by: 
= {B° -V)^ 1 -B°(V -v 1 )-^. (14) 



dt 



The projection of l|14|l onto B governs the evolution of 
the magnetic energy density: 

d B° ■ B 1 v 1 ■ B° |B°| 2 j B ■ B° 

D • V V • V ; . 

An 

(15) 



dt 4tt 



477 47T 



The GW source terms in (|13H15|I are given by: 



B° x d I ~ h * \ . Bid ( h + 

Je = \ h + > Ob = -^ttt: h > 



2 dz 







2 dt 







(16) 



B. Eq. of state & Conservation of number density 

From the first law of thermodynamics, dll — dQ—pdV, 
we can find the internal energy per unit mass, U, as a 
function of the pressure, p, the specific volume per unit 
mass, V = l/p, and the heat flow that, however, vanishes 
under the mhd condition: dQ = 0. We assume an adia- 
batic equation of state p = Kp 1 , where 7 is the adiabatic 
index, which lies in 4/3 < 7 < 5/3. The total relativistic 
matter energy density with respect to the 4- velocity of an 
ideal fluid (with the velocity of light included explicitly) 
is given by: 



/i = p(c 2 + U) = pc 2 + 



P 



7-1' 



(17) 



and the relativistic enthalpy is defined by w° = p° +p a . 
From these expressions we can derive the proper relativis- 
tic sound velocity by considering the change in pressure 



with p at constant entropy: 

, dp 
s dp 



IP 



(18) 



.■ 1 



The matter density p can be solved from the covariant 
conservation of proper number density, n — p/m e , which 
is given by V ' a {nu a ) = 0. In the orthonormal comoving 
frame to first order one finds: 



dt 



(19) 



C. Energy conservation 

The conservation of energy and momentum follows 
from the divergence of the efe (Q) as: 

V b T ah = V 6 • [(p + p)u a u b + pg ab ] - F ab Jb = 0. (20) 

Contracting (|20[1 with U c a leads to the energy conserva- 
tion equation: 



d_ 

dt 



[r 2 (p + p)-p}+W-[T 2 (p+p)v}=0, (21) 



In the comoving frame (r = 1) we can eliminate p in 
favor of p from (|17f) and l|19|) : 



^+ 7 p°VV=0. 
dt 



(22) 



Note that from hereon, (|19fl is redundant and can be 
dropped from our set of equations altogether. It can be 
solved separately to find p 1 from v . 

In the cold non-relativistic limit when the internal en- 
ergy is negligible with respect to the rest-mass energy 
(p -C p), flUl reduces to lfl§|) as in Paper I. 



D. Conservation of momentum 

The equation of motion or momentum conservation in 
a 3 + 1 split is projected out by H c a : 

d 

— [T 2 {p + p)v] + V • [r 2 (^ + p)vv + P I]=jxB, 



(Z+P^+Vp^i'xB (23) 

The Lorentz force couples the matter to the electromag- 
netic fields and the GW source terms, which is apparent 
from combining (|13fl and l|23l) : 



dv 1 , (BO-V)^ 1 [B^-B 1 

B°d{v 1 -B°) j'e x B° 



4vr dt 



4tt ' 



(24) 
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where we have defined: 

IB°I 2 

w tot = w°+ ] —^. (25) 

47T 

This form of the equations of motion is also found 
by explicitly evaluating the divergence of the electro- 
magnetic part of the stress-energy tensor, V&T^ = 
±V b {F a c F bc - \g ab F cd F cd ), instead of using Ampere's 
law. The time projection, V^T 06 , then results in a con- 
servation equation for the total energy density, e.g. i|15|) 
and H22fl combined. 

To summarize, in this section we have obtained a closed 
set of partial differential equations of z, t for the 16 vari- 
ables B, E,j, t, v, fi,p, p that constitute the general rel- 
ativistic mhd description of a GW propagating through a 
magnetized relativistic plasma. We will first solve these 
equations algebraically in Fourier and Laplace space in 
the next section and subsequently derive the space-time 
solutions in Sect. Ivl 



to allow for growing amplitudes as in Paper I. The wave 
equation (|27|l can be written in a symmetric matrix rep- 
resentation since the non-linear general relativistic MHD 
equations form a set of symmetric hyperbolic partial dif- 
ferential equations: 

Dv 1 = j£ w , (28) 

where D is: 

/ u) 2 (l-u 2 A± )-k 2 u\ n -0 2 -fc 2 )u A ||M A± \ 

oj 2 -/c 2 m a|| 

y -(w 2 -ft 2 )l4 A ||«A± W 2 (l-M 2 L|| )-fc 2 (M 2 n -M 2 v|| ) J 

and, for h + x oc e^( z -') as discussed in Sect. Ill Bl (but 
see Sect. I Villi , the GW source terms J GW are: 

ri iuJ 2 u A ± ( " A l^+ \ 

Jgw = — uJ _ k u A\\hx ■ (29) 
w y -UA±h + J 

B. Dispersion relation 



IV. PLASMA WAVES 

In a relativistic, magnetized plasma the proper sound 
velocity of a compressional wave is given in l|18[l and we 
define the Alfven velocity of a non-compressional shear 
wave in the magnetic field, ma, and the velocity of a 
mixed magneto-acoustic wave (msw), u m , by: 



"a 



\B 



(112 



47TW t ot 
+ 



\B 



012 



IP 



wtot 47rw t ot 



(26a) 
(26b) 



To construct a wave equation for the plasma perturba- 
tions we take a second time derivative of H24J) and elim- 
inate B 1 by {HI El and p 1 by {22} • In terms of the 
characteristic velocities (|26al I26bfl we find: 



dt 2 



v 1 - 



d 2 



(v 1 ■ M A ) = 



(m a • Vfv 1 - m a (m a ■ V)V -v 1 + GW terms, (27) 
where the GW source terms are now given by: 



WtQt 

47T 



d 

V (j B • ma) - -t^ (j'e x m a ) - (m a • Vjjl 



A. Symmetric matrix representation 

We are considering GW-propagation along the z-axis 
with the wave vector k = (0, 0, k) at an arbitrary an- 
gle 9 to the ambient magnetic field, which is assumed 
to lie in the x — z plane. To solve the system of differ- 
ential equations algebraically, we Fourier transform with 
respect to time and Laplace transform the spatial part 



A non-trivial solution for the plasma waves in l|28|) 
requires the determinant of D to vanish. Solving for 
k — k z (u), since u> is fixed by the driving GW, we find six 
solutions: 



uj = ±k A u A cos 9 — ±k A u A \\, 



(30a) 



u) = ± 



V2 



Uiul + c 2 ^) yJl±y/fT^j, (30b) 



where we have defined the auxiliary parameter: 



a(0) 



4c s«A|| 



(u^ + c^u 2 ) 2 ' 



(31) 



The negative sign in (|30bjl refers to the the relativis- 
tic proper slow magneto sonic waves with phase velocity 
u s = uj/k s , and the positive sign to the fast magnetosonic 
waves with uf — u>/k{. Together with the Alfven waves, 
we have obtained a 6 x 6 mhd representation that in the 
special-relativistic limit reduces to the equations in [l5j . 

In a low plasma-beta (/3 p i = 4itp/Bq), e.g. strongly 
magnetized, plasma where u A 3> c s , the fast mode re- 
duces to the magneto-acoustic mode f (126 hi) ) slightly al- 
tered by the presence of the gas, whereas for a high 
plasma-beta (c s 3> u A ) it reduces to a sound wave with 
velocity (|18|l . For the slow mode, the behavior is essen- 
tially the other way around. 

The angular dependence is similar: for parallel propa- 
gation (0 = 0), Uf — >u A and u s — >c s , whereas for perpen- 
dicular propagation (9 = n/2) which we studied in Paper 
I, Uf — >u m and it s — >0. 



C. Coupling 

In this section we present a formal derivation of the 
solutions to the inhomogeneous wave equation (|28|l fol- 
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lowing [3. From (EHJ we have: 



v 1 = D- 1 ^ = 



A 



(32) 



where Aij(io, k) and A(oj,/c) are the matrix of cofactors 
of D and its determinant, respectively. These are related 
by: D ki \ kj = SijA. 

Each solution of A(u>, k) = can be identified with a 
wave mode M with w = ui{kM) and w(— fcjw) = — <^(&m)- 
The determinant can be factored into these wave modes: 



E. Slow and fast MSW driven by + polarized GW 

As we expect from our considerations in Paper I, a 
+ polarized GW excites slow and fast magneto-acoustic 
waves in the plasma. The velocity components are: 



vl(k,u) 



2 {to 2 - k 2 u 2 ){uj 2 - k 2 u 2 ) uj-k 



v z (k,uj) 
tan 



k 2 c 2 



, (36a) 



(36b) 



A(w,jfe) = {uj 2 



k 2 u 2 All )(uj 2 - k 2 uf)(cu 2 - k 2 u 2 s ) = 



s«A|| 



i± 



(33) 



The unit polarization vector for a wave in mode M, 
riM(k) with riM • n* M = 1, can be constructed from A^: 



n M i(k)n* Mj (k) 



Xtjjcj, k M ) 
\ii{uj, k M ) ' 



(34) 



D. Alfven waves driven by x polarized GW 

The wave solutions can now be evaluated in Laplace 
space. The \ yy component couples to the x -polarized 
source term and excites Alfven waves, viz perturbations 
of the magnetic field perpendicular to the background: 



i h x Lou^U\ 
~2~ 



uJ - 



k 2 



-vl(k,u>) 



B°. 



LO - 
LO ■ 



ku 2 



(35a) 
(35b) 



and from (|22|l one can easily find the pressure: p 1 (k, u>) — 
~ r yp°v z (k,uj). The magnetic component can be derived 
from (|36|l as: 



v, sin ( 



vz cos t? — 



1-«A v l 



cos 9' 



(37) 



The fact that the magnetic field perturbation is orthogo- 
nal to the direction of propagation of the GW is dictated 
by V • B = 0. 

The polarizations in the magneto-acoustic modes are 
illustrated by Fig. 



* z 




and similarly for E\. z and j\ z . Note that because Figure 2: Orientation of the perturbations in the MSW modes 
the Alfven wave propagates obliquely with respect to 
the background magnetic field, the electric field is no 
longer divergence free and consequently the Alfven wave 

is accompanied by perturbations in the charge density: v. SPACE-TIME SOLUTIONS 

v • /<; ; ,/,•/:.' = 1 ~t ' . ~ ~ 

The polarization of the Alfven wave components is To find the solutions in spaC e-time we apply the inverse 
summarized m Fig. Q] Fourier and Laplace transformations to the results of the 

previous section. We define the phases of the wave modes 
as </>^ = ±k\z — Lot and similarly for ^ , cjr^, and (j) g — 
uj[z — t). 



Alfven waves 




The most straightforward are the Alfven waves i|35|) 
with wavenumber k& — lu/uj^u: 



Figure 1: Orientation of the perturbations in the Alfven mode 
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B2 



4 



-"A|| 



1 + ^A|| 



Ml 



'A|| 



"A|| 



(38) 
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and similarly fon;' (z,t), E* z (z,t),jl z (z,t), andr^z^) 
(see Appendix Since the Alfven waves are non- 

compressional, we do not find a p 1 or \i x contribution. 



B. Magneto-acoustic waves 

Slightly more complicated is the coupled superposition 
of slow and fast wave modes, polarized in the x— z plane. 



h+ u 2 AA 



z 4 u 2 - 

h+ u 2 A± u 2 
4 u 2 — u 2 



l + u s e l fi 
l-u s u s 



l-u s 

1 + U S U s 



4 e #g 



l+me 1 ^ l-u { c l 



1 — Uf Uf 1+Uf Uf 



1 - 

l-u 



and: 



vi tan ( 



4 uf- 



1+ttf e* 



1-ttf e^f 



1 — Uf Uf 1 + Uf Uf 

4^e^ l+use^ l-u s e 1 ^ 7 Au 2 ^^- 



(39) 



(40) 



1 



1-U S U s 



1+U S 



1 



Note that for c s 10, ^TB . If also 9 = ir/2 the limiting 
behavior is: Ut — >ma, u s — >0 and we retrieve our original 
idealized solution. 

The remaining magneto-acoustic wave components, 
B l 0> 0: E l OM), il {z,t), p 1 (z,t), and fi 1 (z,t), are equiv- 
alent superpositions of slow and fast waves but with dif- 
ferent relative amplitudes, summarized in Appendix 1*21 



should also have a transverse component, since the ampli- 
tudes are proportional to B° sin 8 sa B°9 for small angles. 
Explicitly: 

= „ ^cj [e ^] +0[6% (42) 

The excited slow magnetosonic wave is a purely os- 
cillatory wave propagating both in the forward and the 
backward directions. 



VI. A WARM RELATIVISTIC PLASMA WIND 

We now want to consider a warm relativistic plasma 
wind flowing out in the z direction with constant (back- 
ground) velocity /3 and corresponding Lorentz factor 



7tot = 1/ v 1 — P 2 - We are allowed to use simple Lorentz 
transformations in this general relativistic treatment, be- 
cause the whole concept of gravitational waves as small 
perturbations of the background space-time (linearized 
theory) , relies on the fact that we can treat the GWs as a 
field living in flat space-time as long as the scale of cur- 
vature is much larger than the wavelength of the GWs: 
Tl ^S> Agw, which was verified in Sect. Ill Bl 

From now on, all quantities derived in the previous sec- 
tions for the comoving frame will be denoted by primes. 
Since the phase of plane waves, 4> = k a x a , is invariant un- 
der Lorentz transformations, a boost mostly implies addi- 



tion of velocities. We then have 



and 



''Ml 



"A|| ~P 
1-/3uaii 



7 _ l-/3u f 



£_ < _ u 3 -!3 



l-/3u s 



for the phase velocities and: 



C. Growth 

In general, for arbitrary ua,c s , and 9, we always find 
Uf > u s . In the limit where the phase velocity of the 
fast mode approaches the speed of light, U{ | 1, coher- 
ent interaction with the GW is possible. The amplitude 
of the retreating fast wave (oc (1 + Uf)^ 1 ) is negligible 
with respect to the forward wave (oc (1 — Uf) -1 ). The 
forward wave grows linearly with distance because with 
kf = uj/u[ = lj + Ak to first order in Ak we have: 



1 — Uf 



Uf Afc 



HAk z. 



(41) 



This limiting behavior is the same for all magneto- 
acoustic components (see l|A.3jl ). 

For the Alfven waves, the conditions for growing solu- 
tions are less favorable because of the angular depen- 
dence of the resonance condition. To interact coher- 
ently with the GW, the phase speed of the Alfven wave 
(v p h = ua cosO) has to approach the speed of light. On 
the one hand this means that the wave vector of the 
Alfven wave should be almost parallel to the background 
magnetic field, but on the other hand, the magnetic field 



«)' 



T 



7 



(43) 



The boosted background magnetic field is B° = 
{"fB®', 0, B z ') and in the laboratory frame one finds a ze- 
roth order electric field E° = -0 x B° = (0, -pyB%, 0). 



A. Relativistic Alfven waves 



We define 7^ = 1/(1 - u\u) + 0[6 2 } as the Lorentz 
factor associated with the Alfven speed for small 9 and 
boost (|A.1(1 to the laboratory frame as an example: 



4 1-«A||0 



{4[uA,|(l + /3 2 )-/3(l + uin)]< 



1 + «A||\ 2 e^i / 1 + /3 \ 2 e^A 



y 4 



1 + 13 J 7 4 \l + u A \\J 7; 
Or in the ultra-relativistic limit ua T 1 (and (3 ~ 1) 

«$(*,*) «-^«*9[e<*]. 



(44) 



(45) 
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B. Relativistic MSW 

In a Poynting flux dominated force-free plasma wind, 
where the magnetic energy density strongly dominates 
the matter density, the plasma flows out at ultra- 
relativistic velocities. In this regime the phase velocity of 
the fast mode approaches the Alfven velocity, Uf ~ ua T 1 
and the phase velocity of the slow mode becomes negli- 
gible, u s I 0. The solutions are therefore quite similar 
to those found in Paper I, the main difference being the 
angular dependence. For instance: 



we find: 



sin 2 (9 



7 4 7 3 (1 -/3 cos 6) 2 
v' z h + sin 9 (cos 8 ~ (3) 
7^ ~ 4 7 4 (l-/3cos6») 2 



cuz 3 [e^l , (46a) 



luz 3 [e^*] . (46b) 



Appendix HI lists the limiting behavior of all the remain- 
ing MSW components in a relativistic wind. 



VII. DAMPING OF THE GW 

To find an evolution equation for the gravitational 
waves, we project the transverse traceless part of the 
stress-energy perturbation as in J2J (as we discussed at 
the end of Sect. Ill El JIJ and (J2J also hold in the ONf). 

The perturbation of the matter part only has ST^ = 
ST^ — p 1 (z,t) which is purely gauge (the trace can be 
removed by a gauge transformation). Only the magnetic 
field interacts directly with the gravitational waves. All 
the other perturbations are excited through the mhd pro- 
cesses in the plasma. In the general case where oblique 
h+ GW excite slow and fast magneto-acoustic waves and 
h x GW excite Alfven waves, (0 is replaced with: 



(47) 



As we discussed in Sect. Ill Bl we have assumed that the 
GW have a slowly varying amplitude Tt(z). We have ne- 
glected this variation in studying the interaction with the 
plasma since in the short-wavelength approximation all 
the derivatives have It is however this 

slowly varying amplitude that describes the damping of 
the GW. We could find an order of magnitude expression 
for H(z) by integrating gj) using ljA.21) and J^HJ, but 
this wouldn't be entirely self-consistent since in deriving 
(|A.2|I and l|38|) we already assumed that w — k for the 

GW. 

It is, however, possible to derive a dispersion relation 
for the damped GW and the excited MHD waves simulta- 
neously in a self-consistent way. If we assume that the 
GW oscillates at a fixed frequency but leaves the spatial 
dependence unspecified (e.g. a boundary value problem 
with h(z, t) (x h(z)e~ luJt ) we can still solve the full set of 
MHD equations in Laplace space. For the magnetic field 



Bl(k,u>) = 



k 2 U 



A|| 



B° x h+(k,Lu) uj 2 + k 2 u\ 



k 2 u\ ' 



(48a) 



(48b) 



where the second expression is in the limit of a Poynting 
flux dominated wind with ua3>c s J, and h + ^ x (k, lu) are 
the Laplace transforms of h + _ x (z, t) (the reason that lO 
look different from (|35bjl and 13 7|) is that derivatives of 
h-\- t x depend on both u> and k in this more general case). 

If we insert (|48l) in the Laplace transform of Q47p. 
h +tX (k,uj) drop out and we find the self-consistent dis- 
persion relation for the coupled fast magnetosonic - grav- 
itational mode: 



uj 2 - k 2 = 2(B X 



,lu 2 + k 2 u\ 
uj 2 — k 2 u\ ' 



(49) 



and the coupled Alfven - gravitational (hx) mode which 
looks the same but with ua replaced with Ua.\\- 
The two modes allowed by ()49f) are given by: 



k\ 2 u \ 



1 + "A 

l-ul 




m terms of e = — , — % 
(|50|l reduces to: 



(B" 



2u A 



, (50) 



For small e, 



kf = or i 



ok 2 } 



i 2 2 
K 2 U A 



= uj 2 (1 + 2e) + 0[e 2 }. 



(51a) 
(51b) 



Mode 1 corresponds to a superluminal mode with phase 
velocity u>/k\ > 1 and mode 2 is subluminal: 



v. 



Ph,i = 1 + 

,2 

pli-2 _ 1 _ 



= 1 - 



= 1 



(2B» 



with Uph.ifgi.1 = 1 and w p h,2Wg r .2 = u\. Both modes have 
group velocity diojdk\ 2 < 1 as long as (consistent with 

8nG (Bl) 2 



< uj(Ak)c. 



(52) 



Cf (J, 

In the limit e -c 1 , i|51b[l reduces to a fast magnetosonic 
wave in flat space-time (or, equivalently, to the Alfven 
mode for h x ). Mode l|51a|l reduces to a vacuum GW, 
justifying our approximation /i+. x oc e luJ ( z ~ l ) in deriving 
the plasma perturbations in Sect. IIVI 



VIII. INTERPRETATION 

In this section we present an intuitive interpretation of 
our main results: that a + polarized GW excites magneto- 
acoustic waves and a x polarized GW excites Alfven 
waves in a uniform magnetic field. 
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The driving force exerted by a GW on test particles is 
described by ljT21) . These equations can be illustrated by 
force lines in the plane orthogonal to the propagation of 
the GW 8]. Integrating (|12[) twice with respect to time 
results in the well known equations of spatial deviations 
of test masses in an interferometer detector such as LIGO: 



Sx = -(h+x + h x yo), 8y=^-(h x x -h + y ). (53) 

In the MHD limit these test particles are 'glued' to the 
magnetic field lines, so the magnetic field will exhibit the 
same behavior (in fact the presence of the plasma is not 
required, the magnetic field lines can just be viewed as 
parametrized by J53J)). Since the action of the GW is only 
in the x — y plane, and the magnetic field lies in the x — z 
plane we expect that a + polarized GW results in: 

5B x oc^h+B° x , 5B y <x^h+B° = 0, (54) 

and a x polarized GW excites: 

6B X oc \h*B% = 0, 5B v <x~h x B°. (55) 

This is exactly what we found in the mathematical treat- 
ment of the previous section. Perturbations in the other 
directions are higher order effects. 



- Perpendicular B-field 





Figure 3: The MSW (left) and AlfVen mode (right) illustrated 
as an oscillating vector field. The x — y axes in the right figure 
are rotated by n/2 with respect to those in the left figure. 

Fig. gives a schematic illustration of the perturbed 
magnetic field in the two wave modes. The left figure 
shows the vector field (in arbitrary units) for the oblique 
magneto-acoustic wave propagating in the z-direction. 
The perturbations are highly exaggerated to emphasize 
the effect. Amplification of the magnetic field occurs 
when B® is amplified and dilution when B® is suppressed. 
Since B® is constant the total magnetic field has an overal 
wavy pattern. 

Alfven waves are non-compressional and only set up 
a vibration in the field lines perpendicular to the back- 
ground field (By). This is illustrated in the right figure, 
where the axes are rotated to emphasize the y — z plane. 
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Figure 4: Merging neutron star binary 



In a pulsar environment the plasma initially flows out 
along the open field lines but develops into a force-free 
wind outside the light-cylinder in which the toroidal com- 
ponent of the magnetic field dominates the poloidal one. 
Here the magnetic field is predominantly perpendicular 
to the radial propagation of the wind as illustrated by 
Fig. 0J Gravitational waves would mainly excite Alfven 
waves in the former region, whereas in the latter case the 
magneto-acoustic waves are favored. In the relativistic 
wind the GW frequency is red-shifted and the interaction 
is suppressed by 7 -2 , but the interaction length scale 
becomes very large. 



IX. CONCLUSIONS 

We have studied the propagation properties of a plane 
polarized gravitational wave in a magnetized astrophysi- 
cal plasma in the most general case. Both polarizations 
of the GW have been taken into account. Oblique propa- 
gation with respect to the background magnetic field was 
studied including pressure terms, and relativistic veloci- 
ties. The only approximations in this treatment are the 
MHD conditions for the plasma, the linear perturbative 
approach, the geometric optics or small wavelength limit 
for the gravitational waves, and our assumed geometry 
of the interaction region. 

The result is a very rich astrophysical problem, where 
all three fundamental plasma wave modes are excited, 
two of which can interact coherently with the driving 
gravitational waves, and as a result grow linearly with 
distance, dissipating GW energy into the plasma. 

Alfven waves are excited already in linearized theory 
by x polarized gravitation waves propagating at an angle 
with respect to an ambient magnetic field, a result that as 
far as we know has not been found before (note that both 
conditions, oblique GW propagation and x polarization, 
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are required). The Alfven waves are non-compressional 
shear waves that have orthogonal electromagnetic com- 
ponents with a corresponding drift velocity in the plasma, 
a current flowing along the electric field and a deviation 
from charge neutrality caused by the divergence of the 
electric field, but no pressure or density components. To 
interact coherently with a GW one has to find the opti- 
mum of longest interaction length scale and significant 
amplitude as a function of the angle between the Alfven 
wave vector and background magnetic field as was dis- 
cussed in Sect. IV CI 

As we already derived in a Paper I, + polarized grav- 
itational waves excite magneto-acoustic waves propagat- 
ing parallel to the gravitational waves. In this paper we 
have generalized our previous treatment to include the 
oblique magnetic field and a relativistic equation of state 
with non-vanishing pressure. As a result we find both 
the slow and the fast magneto-acoustic waves with phase 
velocities that depend on both the electromagnetic and 
the matter properties of the plasma. 

The magneto- acoustic wave modes are compressional 
waves, that excite pressure, density, and magnetic field 
gradients along their wave vector direction, but no per- 
turbation of charge neutrality. The electric and the mag- 
netic field perturbations and the wave vector are mutu- 
ally orthogonal (E 1 ±B 1 J-k Sy i ), but the drift velocity is 
no longer exactly perpendicular to the magnetic field due 
to the pressure. Of course, the plasma motion in a GW is 
non-compressional but it generates magnetic field com- 
pression if it propagates across a magnetic field either 
in a vacuum or in an ideal frozen-in plasma, and hence 
couples to the magneto-acoustic wave. 

As to the slow MSW, no coherent interaction can occur 
with the GW since the phase velocity of the slow mode is 
always much smaller than that of the gravitational waves. 
Therefore the amplitude does not grow in time or with 
distance and cannot become significant. 

The most effective interaction occurs between the GW 
and the fast magneto-acoustic wave. The phase velocity 
of those waves can easily approach the speed of light in 
a strongly magnetized plasma (in the limit ua 3> c s ) as 
we derived in Paper I and therefore the waves will grow 
linearly with distance. As the plasma waves grow, the 
amplitude of the gravitational waves decays correspond- 
ingly but as long as the total length of the interaction re- 
gion is smaller than the background curvature produced 
by the plasma itself, this will only be a small fraction of 
the total gravitational wave energy. 

Astrophysical applications of the above interactions 
lie in sources of strong gravitational waves embedded 
in strongly magnetized plasmas. Examples are non- 
spherical rotating neutron stars, fast rotating neutron 
stars that are unstable to torsional oscillations, non- 
spherical supernovae collapses, magnetars, and the pro- 
genitors of gamma-ray bursts (non-spherical collapse of a 
massive star or a merging neutron star binary). Most of 
these sources are probably accompanied by an extended 
strongly magnetized and force-free plasma wind, flowing 



out at ultra-relativistic velocities. Since the coupling con- 
stant of gravitational waves is exceedingly small, only in 
the most extreme sources will the interaction with the 
plasma be significant. 

Fig. 0] shows as an example a binary neutron star as a 
gamma-ray burst progenitor. For such a merging binary 
with a magnetar class surface magnetic field of 10 12 T 
[l8|. an angular frequency at the end of the spiral-in 
phase of the order of 1 kHz and a force-free wind flowing 
out with a Lorentz factor of T ~ 100 up to a fraction of 
a light year, a total energy of 10 43 J can be transferred 
from the gravitational waves to the wind. 
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Appendix: ALL SPACE-TIME SOLUTIONS 
1. Comoving Alfven wave components 

The components of the Alfven wave can be derived 
most easily by starting with it's velocity component: 



Vy _ h, 



«A1 



1+UA ll„#t 1_UA II 



4u A ||e I0g 

^17 



(A.l) 



1-uaii 1+«a|| 

From E 1 = —v x B° we find the electric field: 

El{z,t) _ El(z,t) _ v 

whereas the current density follows from Ampere's law: 



1-ui 



A|| 



The charge density is given by V • E 1 = 4-7TT 1 : 
T (z, i) iujh 



Bl tan 9 



■"A|| 



1+wah 



2. Comoving MSW components 

The most complicated MSW component is the magnetic 
field that clearly betrays the mixed nature (gas and elec- 
tromagnetic fields) of this mode. 



Bi(z,t) 



B%, 



4(u 2 -w 2 ) 



[C(z,t) + D(z,t)} 



(A.2) 
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with: 

C(z,t) 



2(u 2 -u 2 )(l + u 2 A )e i ^ Uiiuf+ujje^t 



(l-u 2 )(l-u 2 ) 



1 — Uf 



Ufiuf-ujje^t ^u s (u s + u\)e l,l, ° | u s (u s -u 2 A )e l4 



1+Uf l-u s 



l + u s 



and: 




Note that in the relativistic limit (u s J. 0) D vanishes. 
From £7 = — t; 1 x £?° we have: 

E y (z,t) = B x vl(z,t) - B° z vl(z,t). 

The current density is most easily found as: 



o 



B°cos9 c s 2 
and the pressure is derived from 122J1 



vl{z,t), 



P 



h + u 2 A± u s 



4 u 2 — u 2 



l+w s e l ^ + l-u s e^ 4e^« 



1~~U S U s 1+W S U s 1 — U 2 



h + u AJ _Uf 
4 u 2 — u 2 



l+u f c^f + l-u f e^f 4e^ 



1 — lif Itf 1+Mf Mf 

which readily leads to the energy density: 

1/ ,x ^(M) 

M OM) = — 5 — • 



l-u 2 



3. All ultra-relativistic Alfven wave components 



First we evaluate B' y in terms of laboratory quantities 
(boosted phase velocities etc. j2l| ): 



h -v 2 R° 
47 



■ "A|| \ 2 e^i 



{4[(l + /3 2 )(l+ M i || )-4/3 UA ||]c^ 



1+/3 



1 + 
L + u A || 



e a 

lux 



(A.4) 



Now we can express all other components of the Alfven 
wave in terms of (|A.4J) and l)44[l: 



5 y 


= 7(5; 


+ {3E' X ) 


= J B' y ~ fijSZVy, 






+ PB' y ) 


= 1 pB' y - 7 B x Vy, 




= % = 


id' 1 

~ B x v y z 


= -B X Vy, 




= f x = 


iu> 


B° z Vy 


87r7 2 72 A 


("A 






= Jti'z 


P/?t') = 


B° 

^ X 


lUVy-2jfj—— 

oz 


r 


= 7(r' - 




B° 

^ X 

8tt 


dVy 

iujpvy-2'y—- 
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4. All ultra-relativistic MSW components 



We will only consider the limit of a Poynting flux domi- 
nated ultra-relativistic wind and Lorentz transform i|A.3|) 
(with w = 7(0/ + 13k') ~ 27a;'): 



In the limit of a Poynting flux dominated wind where 
Uf ~ ma T 1 and po I so ii s ~ c s I the fast mode 
can interact coherently with the GW and we find growing 
amplitudes for all components: 

Bl(z,t) vl(z,t) v x (z,t) f{z,t) 



Bq sin ( 

El(z,t) h+ . 

Bn 2 



/i° sin ( 



[ e *«] 



(A.3) 



/' 



7(i?; - /3^) * 2 7 i£ = ^B2w 2 3 [e*«] , 



7 (£;-/3£Q 



—B x , 
sin 2 6> 



yB 8 7 3 5° (1-/3 cos 
, ft 



= 7M = 



/U° sin ( 



47 2 1 - (3 cos ( 



-w 2 z » [e^] , 
9 [e i<fe ] . 
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